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Vibrations  of  Doubly-Rotated-Cut  Quartz  Plates 
With  Monoclinic  Symmetry* 

R.D.  Mindlin 

P.0.  Box  385,  Grantham,  NH  03753 

ABSTRACT:  For  certain  doubly  rotated  cuts  of  quartz,  the  elastic 
stiffness  constants  have  the  same  symnetry  and  absolute  values  as  those 
for  certain  rotated-Y-cuts ;  but  four  of  the  thirteen  constants  have  signs 
reversed.  Mathematical  solutions  of  corresponding  problems  for  the  two 
types  of  cut  have  the  same  form  but  the  numerical  results  may  be  the  same 
or  different  according  as  the  constants  with  changed  signs  enter  the  so¬ 
lution  as  even  or  odd  powers  or  products.  Examples  of  both  are  exhibited. 

I.  Equations  for  Doubly-Rotated-Cut  Quartz  Plates 

Alpha-quartz  has  an  axis  of  three-fold  symnetry,  say  X3  ,  and 
three  axes  of  two-fold  symmetry  one  of  which  is  designated  as  X.j  in  a 
right-handed,  rectangular  coordinate  system  X^  ,  i=l,2,3.  A  doubly 
rotated  set  of  axes  is  obtained  by  rotating  the  X^  system  a  positive 
angle  8  about  X^  and  a  positive  angle  $  about  X,  to  a  new 

investigation  supported  by  the  Office  of  Naval  Research J 
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orlentation  .  The  direction  cosines  of  the  xi  axes  with 

respect  to  the  axes  are 


X1 

x2 

X3 

X1 

£ll 

=  cos  <p 

£^2  =  sin  0 

£13  s  0 

x2 

£21 

-  -sin  0  cos  0 

£22  =  cos  0  cos  6 

£23  =  sin  6 

x3 

£31 

=  sin  0  sin  6 

£32  =  “  cos  4  sin  8 

£33  =  cos  e 

Rotated-Y-cut  and  doubly-rotated-cut  plates  are  cut  with  faces  perpendicular 
to  x 2  >  as  shown  in  Fig.  1. 

The  elastic  stiffness  constants  crs^u»  r,s,t,u  *  1,2,3  (or,  in  the 
reduced  indicial  notation:  cpq,  p,q=1...6),  referred  to  the  rotated  axes 
xi  ,  are  expressed  in  terms  of  the  constants  c?^  ,  referred  to  the  Xi  , 
by 

crstu  =  C1 jkH^ri £s j  Atk£uA  ^ 


simmed  over  i,j,k,£  »  1,2,3. 

Of  the  21  possible  constants  cpq(*  c°p)  »  referred  to  the 
coordinates,  only  six  are  independent  inasmuch  as,  for  a-quartz  [1], 


oooo 
c22  *  C1T  c55  *  c44’ 


0  0  0 
c14  "  c56  *  "c24* 


1 

7 


9 


000000000 
c15  *  c25  *  c35  ■  c45  *  c16  *  c26  *  c36  *  c46  *  c34 


0  . 


(3) 


! 


-3- 


The  values  of  the  remaining  c°  , 

pq 

as  given  by  Bechmann  [2],  are 

=  86.74 

c°2  »  6.98 

c33  "  107‘2 

c°3  =  H.91 

(4) 

cj4  =  57.94 

c?4  =  -17.91 

in  units  of  1010  dyn/crn^  or  10^  N/m2  . 
From  (2)  and  (3)  we  have 


crstu  =  cV£rl£sl£tl£ul+£r2£s2£t2£u2  +  T  (£rl£s2  +  £r2£sl  ^£tl£u2  +  £t2£ul ^ 


*  C33£r3£s3£t3£u3 

+  c44^i!'r2^s3  +  £r3£s2)(£t2: £u3+£t3£u2}  +  (£r3£sl+£rl£s3)(£t3£ul+£tl£u3)] 
+  ci2££r]£sl£t2£u2+£r2£s2£tl£ul  “I  {£rl£s2+£r2£sl)(£tl£u2+£t2£ul^ 

+  c13^£r3£s3(£t1£ul  +  £t2£u2^  +  £t3£u3(£rl£sl+£r2£s2^ 

+  c?4C^r2£s3+£r3£s2^£t1£ur£t2£u2^  +  ^£t2£u3+£t3£u2^£rl£sr£r2£s2^ 

+  (£r1£s2+£r2£sl^£t3£ul+£tl£u3^  +  (£tl£u2+£t2£ul *(£r3£sl+£rl£s3^  * 

(5) 


Finally,  upon  substituting  (1)  in  (5),  we  find 
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cn  *  cn 

c12  “  c°2  Cos2e  +  c13  Sin2e  +  c°4  S1n  20  cos  ^ 

c13  *  c12  s*n2e  +  c13  cos2e  “  c14  sin  20  cos  ^ 

c14  *  ^C13  "  c12^  sin0  COs6  +  C14  C0S  20  C0S  ^ 

c15  ”  c14  cos  0  sin  3 $ 

c16  "  c14  sin  0  sin  3$ 

c22  *  C11  cos4e  +  c33  sin4e  +  (c44  +  7  c?3^  sin2  20  +  4  c°4  sin  0  cos30  cos  3<$ 

^23  *  (cn+c33  “  2  ~  4  c44)  sin  29  +  c^3  +  7  c^4  sin  46  cos  3<f 

c24  =  "  C11  S1n  0  cos3e  +  c33  S1n39  cos  0  +  j  (c44  +  ^-  c°3)sin4e 

-  c°4  cos©  cos  30  cos  3$ 

0 

C25  =  C14^3  Sln  0  “  !)  cos  e  sin  3<$ 

_  A 

C26  “  3  C14S1n  0  c°Se  Sin  3$ 

c33  *  c^  sin40  +  c°3  cos40  +  (cj4  +  ^  c°3)  sin220  +  4  c°4  sin3e  cos0  cos  3$ 
c34  *”  c°i  sin30  cos  0+  c°3  sin  0  cos30 (cj4  +  ^-  c°3)sin  40  -  c°4  sine  sin  39  cos  3 

_  A 

C35  '  "  3  c14  sine  cos  0  sin  3$ 
c36  *  c°4(3  cos2e  -  1)  sine  sin  3<p 

C44  m  T  (c°i  +  c33  ’  2  c°3  "  4  c44^  sin229  +  cj4  +  j  c°4  sin  4e  cos  3* 

_  A 

c45  “  c14^3  cos  0  ’  1)  sine  sin  3$ 

c46  *  ciV3  s1f,2e  *  0  cos  0  sin  3$ 

c55  *  7  *c11  "  c12^  s1r,2e  +  c44  cos2e  “  c°4  sin  20  cos  3<p 

Cgg  *  “  ^  (c^i  *  c 1 2  ”2  c44)  sin  20  +  c^4  cos  29  cos  3$ 

c66  “  7  *c?l  “  c12^  cos2e  +  c44  s1n2e  +  c?4s1n  26  cos  3 ♦ 


(6) 


5- 


These  are  the  Cp^  which  appear  in  the  stress-strain  relations  referred 
to  the  : 


cijk£Skt 


or 


CpqSQ 


(7) 


in  which  the  strains. 


Sij 


or  S  ,  in  terms  of  displacements,  u.. 


are 


S11  =  S1  =  U1 ,1 

S22  =  S2  =  u2,2 
S33  =  S3  =  u3,3 


2S23  =  S4  =  u3,2  +  u2,3 

2S31  =  S5  =  ul,3  +  u3,l 

2S12  =  S6  =  u2,l  +  U1 ,2 


(8) 


Upon  substituting  (8)  in  (7)  and  the  result  in  the  stress-equations 
of  motion: 


Ttj,i  =  p uj  ■  (’> 

we  find  the  displacement-equations  of  motion: 

Dijuj  •  »»i  <>°> 

or 

Dllul  +  °12u2  +  D13u3  *  pul  * 

^21 U1  +  ^22u2  +  ^23u3  s  P ^2  *  (11) 

s  P«  3* 


°3\“l  *  D32u2  *  °33u3 


In  which  the  D^j(=  D^)  are  the  differential  operators 

°11  =  cll3l  +  c6632  +  c5533  +  2  c563233  +  2  c15333l  +  2  C16  3132  * 

°22  *  c2232  +  c4433  +  c663l  +  2  c64333l  +  2  c263l32  +  2  c243233  ‘ 

°33  *  c3333  +  c5531  +  c4432  +  2  c453l32  +  2  c343233  +  2  c35333l  » 

°23  =  c5631  +  C2432  +  c3433  +  ^C23+CW3233 + ^c36+c45^333l  + ^C25+(W31 32 

°31  =  c6432  +  C3533  +  c153l  +  ^C31+C55^333l  +  ^C14+C56*3132  +  ^c36+c45^3233 
°12  =  c4533  +  C163l  +  c2632  +  ^C12+C66^3132 +  ^C25+CW3233  +  ^c14+c56^333l 

(12) 

where  3  .3.  =  S^Sx.Sx.  ,  3?  =  32/3x?  . 

I  J  J  *  ' 

For  traction-free  planes  parallel  to  the  coordinate  planes,  it  is 
required  that 

on  x^  =  constant: 

T11  *  T1  *  C11S1  +  C12S2  +  C13S3  +  C14S4  +  C15S5  +  C16S6  *  0  ‘ 

T12  =  T6  =  C61S1  +  C62S2  +  C63S3  +  C64S4  +  C65S5  +  C66S6  *  0  ’  (13) 

T13  *  T5  *  C51S1  +  C52S2  +  C53S3  +  C54S4  +  C55S5  +  C56S6  =  0  ; 

on  x2  =  constant: 

T21  =  T6  =  C61S1  +  C62S2  +  C63S3  +  C64S4  +  C65S5  +  C66S6  =  0  * 
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on  x3  =  constant: 


T31  =  T5  =  C51S1  +  C52S2  +  C53S3  +  C54S5  +  C55S5  +  C56S6 
T32  “  T4  =  C41S1  +  C42S2  +  C43S3  +  C44S4  +  C45S5  +  C46S6 
T33  =  T3  =  C31S1  +  C32S2  +  C33S3  +  C34S4  +  C35S5  +  C36S6 


0  , 
0  , 
0  . 


(15) 


II.  Rotated-Y-Cuts  vs_  60°  Doubly-Rotated-Cuts 

CASE  A.  If  $  =  0  ,  8  3*0  (the  rotated-Y-cuts)  then  sin  3<)>  =  0  ,  cos  3$  =1 

and,  from  (6), 


C15  =  C16  =  C25  =  C26  =  C35  =  C36  =  C45  =  C46  =  °‘ 


(16) 


The  remaining  13  constants  are  those  for  monoclinic  symmetry  with  x-j  the 
digonal  axis. 

CASE  B.  If  <|>  =  60o,  6^0  (doubly-rotated  cuts),  then  sin  3$  =  0, 
cos  3$  = -1  and  (16)  again  hold  so  that  the  symmetry  is  the  same  as  for 
rotated-Y-cuts.  Even  if  0  is  the  same  in  A  and  B,  all  the  surviving 
constants  in  B  (except  for  c^  which  remains  fixed)  are  different  from 
the  corresponding  ones  in  A  as  the  last  term  in  each  cpq  has  its  sign 
reversed.  However,  if  0  in  B  is  the  negative  of  6  in  A  ,  nine  of 
the  constants  are  the  same  for  the  two  cuts  and  the  remaining  four  have 
the  same  absolute  values  in  A  and  B  but  are  of  opposite  sign. 


T 
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To  sumarize  the  properties  of  the  two  sets  of  constants 


'pq 


($  »8)  I 


cpq(°»8)  =  cpq{60° ,-0)  for  pq  =  11 ,12,13,22,23,33,44,55,66, 
cpq(°,e)  =  -cpq(60°,-e)  for  pq  =  14,24,34,56  . 


The  displacement  equations  of  motion  reduce  to 


^C1131 +c6632  +  c5533  +  2c56?233^ul  +  ^c12+tW3 132  +  ^c14+c56^3  33  1  ^u2 

+  C(ci3+c55)333i  +  ^c]4+c56^3 l32^u3  =  p U1  * 

^c12+cW3132  +  ^C14+C56^3331^U1  +  ^C2232  +  C4433  +  C663l +2  c243233^u2 
+  [c563-,  +  c2432  +  c3433  +  (c23+c44)3233]u3  =  pu2  , 

^c12+(W3132  +  ^c14+c56^3331^u1  +  ^c563l  +  c2432  +  c3433  +  ^c23+c44^3233^u2 
+ ^c3333  +  c5531  +  c4432 +2  c343233^u3  =  p u3  ’ 

(18) 

and  the  conditions  (13),  (14),  (15)  for  traction-free  boundaries  reduce  to: 


on  x-j  =  constant: 


T11  =  T1  =  c11ul,l  +  C12U2,2  +  c13u3,3  +  c14^u3,2  +  u2,3}  =  °» 

T12  =  T6  =  c56^ul ,3  +  u3,l*  +  c66(u2,l  +  ul,2^  =  °» 

T13  =  T5  =  c55^ul ,3  +  u3,l  ^  +  c56^u2,l  +  Ul,2^  =  0; 


(19) 


on  X£  =  constant: 


T21  =  T6  =  c56(ul,3  +  u3,l^  +  c66(u2,l  +  ul,2^  =  °» 

T22  =  T2  =  c12ul ,1  +  c22u2,2  +  c23u3,3  +  c24(u3,2  +  u2,3*  = 

T23  =  T4  =  c14ul  ,1  +  c24u2,2  +  c34u3,3  +  c44(u3,2  +  u2,3}  =0; 


(20) 


on  Xj  =  constant: 

T31  =  T5  =  c55(ul,3  +  u3,l}  +  c56(u2,l  +  ul,2>  =  °* 

T32  =  T4  =  c14ul  ,1  +  c24u2,2  +  c34u3,3  +  c44(u3,2  +  u2,3)  =  0’  (21) 

T33  =  T3  =  c13ul,l  +  c23u2,2  +  c33u3,3  +  c34(u3,2  +  u2,3*  = °* 

Any  solution  of  the  equations  of  motion  (and  boundary  conditions,  if 
any)  referred  to  rotated  axes  with  $  a  0,  9  =  0'  is  the  same  solution,  at 
least  in  form,  referred  to  axes  with  $  =  60°,  9  =  -0'  .  Whether  or  not 

the  solutions  are  the  same  numerically  depends  on  the  occurrence  of  c^, 
c24*  c34’  c56  as  even  or  oc^  P°wers  or  Products  in  the  resulting  formulas. 
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In  the  following  Sections,  we  review  solutions  obtained  previously 
for  Cpq(0,e)  and  determine  if  the  transition  to  cpq(60°,  -6)  changes 
the  numerical  results. 


III.  Plane  Waves  in  a  Plate 

In  a  plate  with  faces  at  x2  =  +  b*  we  consicier  waves  propagating 
in  the  direction  of  the  two-fold  axis  of  symmetry  x-j  : 

u-j  -  A]  sin(S2x2  +  £3x3)  sin(?1x1  -  ut 

u2  =  -  A2  cos ( ^2*2  +  ^3X3^  C0S(^1X1  -  '  (22) 

u3  *  -  A3  cos(£2x2  +  €3^3)  cos(C-|X1  -  uh.)  . 

Upon  substituting  (22)  in  (19)  and  setting  the  determinant  of  the  coeffi¬ 
cients  of  the  Aj  equal  to  zero,  we  find  the  equation 

-  «</i  -  «•  xi3-  <23> 

in  which  6 . .  is  the  Kronecker  delta, 

*  J 

*11  *  ?ll+e  +  2  ^560r  ’  *23  *  ?56+?24e2  +  ^34^  +  ^23+^44^6r’ 

*22  “  1  +^226  +^44^2  +  2  ^24Br*  *31  *  ^14+^56^  +  ^13+^55^r  ’  (24) 

*33  “  ^55+^44B  +^33^+ 2  ^34er’  *12  s  +  ^l4+*Wr’ 


f 
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Cpq  =  Cpq/C66  * 


6  =  C2/C-, 


2  2-2  2  2 

r  =  ,  v  =  nvq  =  /c66q 


(25) 


In  (25),  6  and  r  are  the  ratios  of  the  wave  length  along  to  the 

wave  lengths  along  x ^  and  x3  »  respecti vely;  V  is  the  ratio  of  the 
velocity  to  the  velocity  v  =  (cgg/p)*;  C-j  (=  2^-j b/rr)  is  the  ratio  of  the 
thickness,  2b,  of  the  plate  to  the  half-wave-length  along  x-j ;  and  n 
is  the  ratio  of  the  circular  frequency  w  to  the  frequency  irv/2b  . 

For  given  6  and  r  ,  (23)  is  a  bicubic  in  the  velocity  ratio  V  : 


V6  +  BV4  +  CV2  +  D  =  0  , 


(26) 


in  which 


B  = 

-  (An  + 

A22  +  A33^  * 

C  = 

A22A33  + 

A33A11  +  A1 1 A22 

0  = 

X  A2  ♦ 
A11A23 

2  2 
A22A31  +  A33A12 

A11A22A33  ”  2  A23A31A12  * 


|  The  coefficients  of  the  bicubic  are  different  for  c-i4»c24»c34*c56 

!  positive  and  negative.  Hence,  for  given  0  and  r  ,  the  roots  of  (26) 

yield  different  sets  of  velocity  ratios  'for  the  rotated-Y-cut 

with  Cpq(0,e)  and  the  doubly-rotated-cut  with  cpq(6O°,-0)  .  An  ex¬ 
ample  is  illustrated  in  Fig.  2  in  which  either  0  is  the  abiscissa  and 

i 

r=  10  or  vice  versa.  In  either  case,  the  lowest  velocity  ratios  V-, 

I  J 

exhibit  little  difference  for  the  two  cuts  —  and  this  is  the  branch  which 

j 

l 

i 


1 

/ 
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would  contribute  predominantly  to  the  fundamental  thickness-shear  mode  of 
the  plate.  However,  the  differences  for  the  upper  velocities  are  larger, 
at  least  for  large  0  and  r  —as  much  as  13%  for  0  =  r  =  10.  These 
differences  survive  any  boundary  conditions  that  may  be  applied. 


IV.  Ekstein's  Solution 


It  will  be  observed,  in  Fig.  2,  that  the  velocity  ratios  are  the 
same  for  Case  A  and  Case  B  if  r  {or  0 )  is  zero.  This  is  the  situation 
for  modes  with  straight  crests  along  x^  (or  X2).  In  the  case  F  =  0  , 

X23  and  change  sign,  in  the  passage  from  case  A  to  Case  B  ,  but 
they  enter  the  coefficients  of  the  bicubic  (26)  only  as  their  product  and 
as  squares  --  resulting  in  no  change  in  roots.  To  examine  whether  this  per¬ 
sists  after  the  introduction  of  free  faces  of  the  plate,  we  consider  Ekstein's 
solution  [3]  for  modes  with  straight  crests  along  x3  in  a  plate  with  free 
faces  on  x2  =  . 

With  T  =  0  ,  and  fixed  and  V,  (23) yields  three  roots 

p 

0n  ,  n=l,2,3.  Thus,  for  steady  state  vibrations,  (22)  may  be  written  as 


u2 

u3 


Aln  sin  5l6nx2  sin  ^1*1 
A2„  cos  5,B„x2  cos  E,X, 

A3n  cos  W'2  cos  51X1 


eiut 


icot 

e 


(28) 


'  — s 
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Then  the  boundary  condi  .''ns  (20): 

T2j  =  °*  J=1*2»3»  on  x2  =  i  b  »  (29) 


result  in  Ekstein's  frequency  equation;  which  may  be  written  in  the  form  [4] 


where 


luin!  =  0,  i  ,n=l  ,2,3, 


uln 

^2n 

p3n 


C12Lln 

c14Lln 


+  L2n  +  256L3n)  cot  Bn5lb  ’ 
+  Bn(c22L2n  +  e24L3n*  ’ 

+  Bn(c24L2n  +  £44L3n)  » 


(30) 


(31) 


Lin  '  cof<Xn1  -  5niVn>''cof(1nn  -  Vn>  •  (32> 

In  the  passage  from  cpp(0,6)  to  cpp(60°,  -0),  the  lin  and, 
hence,  the  (which  depend  on  ci 4»c24 *c56^  chan9e  S19n  for  subscripts 

13,31,23,32  while  the  remaining  terms  in  (31)  and  (32)  do  not  change.  But 

those  uin  which  do  change  appear  only  as  product  pairs  in  (30)  and,  hence, 
the  roots  of  (30)  do  not  change.  These  roots  are  usually  depicted  graphic¬ 
ally  as  a  many  branched  dispersion  relation  between  n  (as  ordinate)  and 
^  (as  absicissa): 

0  =  n(^)  (33) 


I 


-14- 

as  illustrated  in  [4].  Alternatively,  the  abscissa  could  be  1/^  : 

«  =  fld/C,).  (34) 

Suppose  the  plate  has  additional  bounding  planes  Xi  *+  a  at  which  the 
conditions  are  uniformly  point-mixed,  e.g.  vanishing  U£,  T11’T1  2  correspond¬ 
ing  to  "simply  supported"  in  the  elementary  theory  of  flexural  vibrations  of 
plates.  For  real  roots  of  (30),  these  conditions  are  satisfied  by  ^=rmr/2a, 
where  m  is  an  even  integer;  so  that,  for  real  roots,  the  dispersion  rela¬ 
tion  converts  to 


fi  =  fi(a/mb)  .  (35) 

Elimination  of  m  from  the  abiscissa  requires  only  that  each  branch  of  the 
dispersion  relation  (35)  be  replaced  by  a  sequence  of  branches  obtained  by 
multiplication  of  its  absicssa  by  a  sequence  of  integers.  In  this  way,  the 
branches  of  the  dispersion  relation  for  the  infinite  plate  are  converted  to 
the  branches  of  the  frequency  spectrum,  R  vs  a/b,  of  the  "simply  supported" 
plate.  As  the  process  does  not  involve  c^.Cg^  and  Cgg  anew,  the  fre¬ 
quency  spectrum  is  not  altered  by  a  change  of  cpp(0,6)  to  cpC)(60o,-e)  • 
There  is  no  closed  solution  of  the  three-dimensional  equations  for 
the  case  of  free  boundaries  at  x^j  *  + a  and  the  situation  there  is  not 
obvious  inasmuch  as  c^  and  c^g  enter  into  the  traction-free  conditions 

T11  8  T12  ”  T13  8  0  on  xl“ia  (36) 


as  may  be  seen  in  (19). 


V.  Effect  of  Free  Edges 


As  a  substitute  for  the  unavailable  extension  of  Ekstein's  solution 
of  the  three-dimensional  equations  to  accommodate  a  pair  of  parallel,  free 
edges,  there  exists  a  solution  of  two-dimensional  approximate  equations  [5]. 
For  the  case  of  straight  crested  flexural  waves  travelling  in  the  direction 
of  x-j  in  a  plate  with  free  faces  at  x2  =  +  b,  the  three  dimensional 
displacements  are  approximated  by 


u-j  =  x2;1  (x1  )ei")t  ,  u2  =  cJ2(x1)eiu;t  , 

"3  ‘ 

(37) 

and  the  differential  equations  governing  them  are 

K  C56U3,11  +  k2c66(U2,11  +  ^1,1 }  * 

2„ 

-  pw  U2  , 

C55U3,11  +  KC56(U2,11  +  *1,1*  = 

2„ 

-  pu  Uj  , 

(38) 

Y1 1^1 ,11  ”  3b  tKC56U3,l  +i<2c66^U2,1  ^1  ^  = 

2, 

•  1 

where 

<  -  n  /12 ,  ^n  =  cH”cl2^c22”^c14”c12c24^^c44”  c24^c22  ^  *  (39) 


There  is  no  change  of  sign  of  with  change  of  sign  of  c-j^  and  c24; 
so  only  c,jg,  in  (38),  changes  sign  with  the  passage  from  cpc?(0,e)  to 
cp„<60°,-e) . 
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The  displacements  are  taken  as 


U2  *  A^b  sin  Sx-j  ,  U3  =  A^b  sin  Cx-|  ,  =  A4  cos  £xi  •  (40) 

Then,  from  (38), 


(?  -  3si2)a2  +  e56£2A3  +  |a4  =  0  , 

c56C2A2  +  (c5552  -  3n2)A3  +  c56tA4  =  0  ,  (41) 

U2  +  S6^A3  +  (^11^2  +  1  "  “2)A4  =  0  » 


where 


t  =  0>,  c55  =  c55/tc2c66  ’  c56 

The  determinant  of  the  coefficients  of  the 
is  the  equation 


c56/kC66  ’  Y11  =  Yll/3<  c66  * 

Ai  in  (41),  set  equal  to  zero. 


A  /A  A?  \  Afi  2  r  <\A  /.  *  »  a  A>2  -|^4 

Y11  c55  "  c56^  "  n  [3-ri^1  +c55)  +  c55  -  c56]? 


+  3fl2[n2-c55(l -n2)  +  3?11n2  +  c2K3C2  +  9«4(1-fi2)  s  0,  (43) 


rlV 


'56- 


A>2 

which,  for  a  fixed  frequency  ratio  n  ,  is  a  bicubic  in  £  whose  roots 
are  Independent  of  change  of  sign  of  c-j 4 *c24 ,c56  *  ^hus»  as  *n  the  three_ 
dimensional  case,  the  dispersion  relation  does  not  change  with  passage  from 
cpq(0,e)  to  cpq(60°,-e)  . 


1 
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•y 

For  each  fl  ,  (43)  has  three  roots  l*  ,  n=l,2,3,  and  (41)  has 
three  sets  of  amplitude  ratios  A^rA^rA^  .  Let  Sn  ,  n=l,2,3,  be  the 
value  of  A4  for  the  ntn  root  Z*  ;  and  let 


-  r-  -  +  3n2-c55^)/An  , 


2n  "A.  ^n'  56^n 


*3n  *  ^  =  3  ^  ^  ’ 

An  *  (^-382)(c55?n2.3«2)-iL4 


"56^n 


(44) 


for  each  root  .  Then  (40)  may  be  written  as 


u3  =  b 


■3 

*l>°2»  sf"  Vl  • 

3 

!y  Van  sin  Vi  • 


■3 


COS  £nx-j 


(45) 


The  conditions  for  free  edges  at  Xi  *+  a  are:  the  horizontal  and 
vertical  shears,  Ng  and  Q-|  ,  and  the  bending  moment,  M^,  vanish.  Thus, 


on  *  + a  , 
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N5  =  2b^c55U3,1  +  K  CWU2,1  +  *1^  =  0  * 

Ql  *  2bic  tc56U3,1 +kC66^U2,1  +  ^1^  =  °*  (46) 

M,  =  (2b3/3)Yll^ul  =  0. 

Upon  substituting  (45)  into  (46),  we  obtain 

3 

^  fln5ln  cos  5na  =  0  ’ 

3 

5-J  5n“2n  c°s  5„a  *  0  • 

3 

y.  S  e  sin  e  a  =  0  , 
n^t  n  n  n 

where 

aln  =  c55a3n^n  +  K  c56(a2n^n  +  ’ 
a2n  =  c56a3n^n  +  K  c66^a2 n^n  +  1 *  • 

The  frequency  equation  is  obtained  by  setting  the  determinant  of  the 
coefficients  of  the  An  in  (47)  equal  to  zero: 


A-j  tan  ^  a  +  Ag  tan  ^  +  A3  tan  £ja  *  0 


(49) 
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where 


*1  “  Ma12a23  "  a22a13*  ’ 
A2  =  ^2^a13a21  "  a23all^  * 

A3  =  Ma11a22  "  a21a12^  ‘ 


(50) 


Upon  substituting  (44)  in  (48)  and  the  result  in  (50),  the  frequency  equation 
(49)  becomes 


tan  ^a1  ^(£3-^)  tan  +  Z3&3(tf  - 1\)  tan  ^  ■  0  :  (51) 


an  equation  which  does  not  change  when  cpp(0,e)  is  replaced  by  cpq(60°,-e). 


VI.  Vibrations  of  a  Strip 

An  exact  solution  of  the  three-dimensional  equations  exists  for 
coupled  thickness-twist  and  face-shear  modes  of  vibration  in  a  rotated-Y-cut 
strip  with  a  parallogrammic  cross-section  and  all  four  faces  free  of  traction 
[6],  The  displacements  are  u2 *  u3 *  0  and,  omitting  a  factor  e1wt  , 


u1  *  A  sin  £2x2  cos  ^3^56x2  ”  x3)  +  B  S1n  ?2X2  sin  ^3^c56x2_x3^ 

+  C  cos  $2x2  cos  £3(cggX2  -  x3)  +  D  cos  £2x2  sin  C3(c5gX2-x3)  , 

(52) 

where  c5g*c5g/Cgg,  as  before  in  (25). 


The  equations  of  motion  (18)  are  satisfied  if 


Plj2  =  c6652  +  Y55?3  *  y55  =  c55  “  c56/c66  *  ^ 

and  the  faces  at  X2*  +  b  satisfy  the  traction-free  conditions  (20)  if 

2^2^  *  nm  (54) 

where  m  is  an  odd  integer  for  solutions  A  and  B  and  an  even  integer  for 
solutions  C  and  0  . 

A  pair  of  planes  parallel  to  the  x-j-axis,  making  dihedral  angles  a 
with  the  x^  -  x2  plane  and  distant  2c  cos  a  apart,  as  illustrated  in  Fig.  3, 
are  free  of  traction  if 


and 


o  *  arc  tan  Cgg 
2£3c  *  nir 


(55) 

(56) 


where  n  is  an  even  integer  for  solutions  A  and  C  and  an  odd  integer  for 
solutions  B  and  D  . 

The  frequencies  are 
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When  cpq(0,e)  changes  to  cpq{6O°,-0)  ,  the  frequencies  do  not 

o 

change  as  c^g  enters  as  c^g  ;  but  the  mode-shape  (52)  changes  and  a  , 
In  (55),  is  reversed  in  sign  so  that  the  cross-section  changes,  as  illus¬ 
trated  in  Fig.  3.  The  values  of  +  a  for  the  full  range  of  values  of  9 
are  illustrated  in  Fig.  4. 


A 
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Captions  for  Figures 


Figure  1:  Rotated-Y-cut  and  doubly- rotated-cut  quartz  plates. 

and  Xj  are  digonal  and  trigonal  axes  of  symmetry, 
respectively. 

Figure  2:  Comparison  of  wave  velocities  for  Cpq(O»0)  and  cpq(60°,-e) 

as  functions  of  the  ratios,  6  and  r  ,  of  wave  lengths  in 
the  X£  and  x^  directions  to  the  wave  length  in  the  direction 
x^  of  the  wave  normal. 

Figure  3:  Cross  sections  of  strips. 

Figure  4:  Variation  of  dihedral  angles,  a  ,  between  face  and  edge 

planes  of  strip  for  cpC}(0 ,6)  and  cp£}(60°,-e)  as  functions 
of  6  . 
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